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Abstract 

Optical hyperthermia systems based on the laser irradiation of gold nanorods seem to be a promising tool in the 
development of therapies against cancer. After a proof of concept in which the authors demonstrated the efficiency of 
this kind of systems, a modeling process based on an equivalent thermal-electric circuit has been carried out to determine 
the thermal parameters of the system and an energy balance obtained from the time-dependent heating and cooling 
temperature curves of the irradiated samples in order to obtain the photothermal transduction efficiency. By knowing this 
parameter, it is possible to increase the effectiveness of the treatments, thanks to the possibility of predicting the response 
of the device depending on the working configuration. As an example, the thermal behavior of two different kinds of 
nanoparticles is compared. The results show that, under identical conditions, the use of PEGylated gold nanorods allows 
for a more efficient heating compared with bare nanorods, and therefore, it results in a more effective therapy. 

Keywords: Gold nanoparticles; Laser irradiation; Optical hyperthermia; Photothermal transduction efficiency; Temperature 
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Background 

Nanoscale materials have been broadly studied in recent 
years, thanks to their unique optical properties and their 
great potential in the development of biomedical applica- 
tions. One of the most interesting areas is the use of 
plasmonic nanoparticles to enhance the diagnostic and 
treatment methods available for cancer. In this field, au- 
thors such as Letfullin and co-workers have recently de- 
scribed the optical properties, the kinetics of heating 
and cooling, and the spatial distribution of temperature 
of this kind of nanoparticles, providing a better under- 
standing of these processes [1-3]. 

Specifically, rod-shaped gold nanoparticles, known as 
gold nanorods (GNRs), have a great potential because of 
their large light absorption and scattering cross sections in 
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the near-infrared (NIR) region [4]. Thanks to this optical 
behavior, GNRs are able to transform the absorbed energy 
into localized heat. This optical effect is used to develop 
cancer therapies as photothermal tumor destruction either 
by direct enough increase of temperature or indirectly by 
co-adjuvant drugs, at the same time delivered by the par- 
ticle, or already present and activated by the heating [5-8]. 

Our research group has recently developed an optical 
hyperthermia device based on irradiation of GNRs with a 
continuous wave (CW) laser in order to induce in vitro 
death of human brain astrocytoma cells (1321 Nl) [9]. 
Unlike many high-energy pulsed lasers that generally lead 
to particle structure changes and ablation in a very short 
time, CW lasers allow heat dissipation from particles to 
surrounding medium (via phonon-phonon relaxation), so 
they are an appropriate choice in order to use the pro- 
duced heat for the cure of cancer [10]. 

The effectiveness of the developed method was deter- 
mined by measuring changes in cell viability after laser ir- 
radiation of cells in the presence of GNRs. In accordance 
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to Other results in comparable experiments [11-13], ours 
indicated that continuous laser irradiation in the presence 
of the particles induced a significant decrease in cell via- 
bility, while no decrease in cell viability was observed with 
laser irradiation or incubation with GNRs alone. Due to 
the limited capacity of laser penetration in tissues, this 
method could be used in clinical practice as an additional 
aid to surgery for removing brain tumors completely. 

After this proof of concept, our objective was focused 
in getting a better understanding about the working 
principles and physical behavior of optical hyperthermia 
devices. It is not very common to find studies including 
a comprehensive characterization about the global phe- 
nomena in optical hyperthermia systems. Moreover, al- 
though now there are a huge variety of noble metal 
nanoparticles that can be used to carry out this kind of 
therapy, an absolute control about their behavior still 
does not exist. Therefore, it is necessary to develop a 
series of characterization and modeling processes to in- 
crease the effectiveness of the hyperthermia treatments, 
thanks to the prediction of the system response. 

With this aim, a method to calculate the thermal param- 
eters of the system and the photothermal transduction ef- 
ficiency for different kinds of nanoparticles has been 
developed. This method, which allows an easy and effect- 
ive thermal characterization and so predicts the thermal 
behavior of the system, is not only valid for our device but 
also for any kind of optical hyperthermia system. 

Methods 

Optical hyperthermia device 

The CW laser (MDL H808, PSU-H-LED power source, 
Changchun New Industries, Changchun, China) works 
at 808 nm, with a maximum output power of 5 W, 
beam height from base of 29 nm, beam diameter at 
aperture of 5 to 8 mm, and dimensions of laser head of 
155 X 77 X 60 mm^. 

We have used two different kinds of commercial 
GNRs in order to compare their photothermal trans- 
duction efficiency. Both are tuned to the laser source 
and have their maximum surface plasmon resonance 
(SPR) at 808 nm (longitudinal band). 

The first commercial GNRs used are bare GNRs (B- 
GNRs) A12-10-808-100 Nanorodz (Nanopartz, Salt 
Lake City, UT, USA). B-GNRs are dispersed in deion- 
ized water (DI-H2O) with <0.1% ascorbic acid and 
<0.1% cetyltrimethylammonium bromide (CTAB) sur- 
factant capping agent. B-GNRs have an axial diameter 
of 10 nm and a length of 41 nm. 

The other commercial GNRs used are PEGylated 
GNRs (PEG-GNRs) PEG-10-808-50 (Nanoseedz, China). 
PEG-GNRs are functionalized by thiol-terminated meth- 
oxypoly(ethylene glycol) (mPEG-PH) and are also dis- 
persed in DI-H2O. The dimensions of PEG-GNRs are 



equal to the dimensions of B-GNRs (axial diameter = 
10 nm, length = 41 nm). 

The laser is connected to the system via a multimode 
optical fiber with a core diameter of 600 [im, a length 
of 1.5 m, and a power transmission of 90% to 99% 
(600-nm MM fiber, Changchun New Industries, China). 
The laser light from the fiber irradiates the samples 
through a collimation lens (78382, Newport Corpor- 
ation, Irvine, CA, USA), which is in direct contact with 
a 4-well plate containing the samples, which have a 
total volume of 500 |il, and is located on a Teflon sup- 
port. A temperature sensor (FlOO Precision Thermom- 
eter, Automatic Systems Laboratories, Redhill, UK) is 
fixed vertically with the aid of a tripod stand and a bur- 
ette clamp and remains in contact with the samples 
during the experiments (Figure 1). 

Thermal parameters 

In order to determine the parameters that characterize 
and describe the thermal behavior of our hyperthermia 
device, it is needed to develop a thermal model, which 
can be raised from the resolution of an equivalent elec- 
tric circuit (Figure 2). 

In this circuit, P is the delivered power, T{t) is the 
sample temperature which is time dependent, and Cd 
(W/K) and Ct (J/K) are the thermal conductance and 
the thermal capacitance of our experimental enclosure, 
respectively. 

Solving the circuit, we can formulate the equation 
that describes the power distribution, obtaining that the 
delivered power (P) is equal to the sum of the stored 
power in the capacitor (Ps) and the dissipated power in 
the resistor {P^): 

P = Ps+Pd 

= c /(^(^)-f-f-»') + (r(^)-r,ef-„)c, (1) 

at 

In this equation, T^f - m is the reference temperature (the 
subscript m refers to the thermal model), that is to say, the 
initial temperature of our sample before the laser irradi- 
ation that should match the environment temperature. 
Once the capacitor is charged, the temperature had reached 
their maximum value {T^^ _ ,„), also called stabilization 
temperature, which remains constant until the end of the 
irradiation. Under this condition, we can rewrite the 
previous equation removing the time dependence as 

P — {P max-m ~ Pref-m ) C^jf — A P m Cd (2) 

In a conventional experiment of optical hyperthermia, 
a laser source irradiates a sample containing a colloidal 
suspension of GNRs which act as little heat sources. In 
the proposed model, the GNRs are replaced by an elec- 
tric resistor {R) which is connected to a voltage source (V) 
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SO that we dispose of a single heat source delivering a 
known power: 



1^ 



(3) 



The resistor heats up the sample until the stabilization 
temperature (Tmax - m) is reached, and then, the voltage 
sample is shut off and the resistor is immediately removed 
from the sample in order to obtain the cooling curves 



T(t) 




Ct 



Cd 



Figure 2 Electrical equivalent circuit used to obtain the thermal 
parameters of the optical hyperthermia device. 



(which correspond to the discharge curves of the capaci- 
tor) that characterize our experimental enclosure without 
the influence of the resistor that is still kept warm. 

By adjusting these cooling curves to the corresponding 
decreasing exponential equation, we can obtain the 
cooling time constant, which depends on the thermal 
capacitance and the thermal conductance of our system: 



t{s) 



Ct (J/K) 
Cd (W/K) 



(4) 



Thus, from a known power and from the values of r 
and lS,Tm experimentally obtained, we can calculate the 
thermal parameters (i.e., thermal capacitance and ther- 
mal conductance) that characterize the sample enclosure 
of the described optical hyperthermia device. 

In this case, we chose a resistor of 15.2 Ci, the voltage 
source values were 1.5, 2.0, and 2.5 V, and the tested 
sample volumes were 500, 750, and 1,000 \A. We obtained 
three heating and cooling temperature curves for each 
possible configuration. 

Photothermal transduction efficiency 

From the Mie theory and taking into account different 
parameters such as the nanoparticle size and shape, the 
refractive index of the surrounding medium, and the laser 
wavelength, authors such as Zharov describe the optimal 
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Figure 3 Temperature curve of heating, stabilization, and 
cooling, obtained from the thermal model and the 
proposed procedure. 



conditions for nanoparticles to obtain effective laser 
heating in optical hyperthermia applications [14]. 

On the other hand, we can find in the literature ad- 
vanced models that completely describe the heat transfer 
behavior from the surface of nanoparticles presenting the 
heat sources produced by nanoparticles in the spherical 
volume of biological tissue [15,16]. These methods allow 
for predicting the complete thermal response for applica- 
tions to future cancer therapies as nanophotothermolysis 
and nanophotohyperthermia, but we propose a simpler 
approach in order to rapidly compare the photothermal 
response of nanoparticles in optical hyperthermia de- 
vices to be able to select those nanoparticles that allow 
us to obtain better results in each planned therapy. 

With the aim of obtaining the photothermal trans- 
duction efficiency of the two different kinds of GNRs 
used as heat sources in the described optical hyperther- 
mia device, it is needed to measure the time-dependent 
heating and cooling temperature curves of the irradi- 
ated samples containing the suitable colloidal disper- 
sion of each kind of GNRs. Roper et al. [17] determine 
the energy balance used to describe this process 
(Equation 5): 



J2 '"'C/ 



dT 



Qi+Qo-Qe 



(5) 



In the previous expression (Equation 5), m and Cp are 
the mass and the heat capacitance of each component of 
the irradiated sample, respectively, T is the temperature 
of the sample, Qj is the calorific energy that GNRs gen- 
erate (energy source), Qo is the baseline energy of the 
sample (represents the temperature rise of the sample 
due to the direct heating of the laser source), and Qext 
represents the energy flux transmitted out of the irradi- 
ated area. 

The term Qj represents the heat that is generated due 
to the electron-phonon relaxation of plasmons in the 
surface of GNRs that takes place because of the irradi- 
ation of the particles at the SPR wavelength A: 




0123456789 10 

Figure 4 Relation between P (W) and AT^ (K). Lineal fits for eacli 
tested value: 500 |jl (blue), 750 |jl (red), and 1,000 |jl (green), whose 
values of thermal conductance are 0.052, 0.052, and 0.048 W/K, 
respectively, fl^ is the average squared error of each fit. 



Qi = 7(1-10- 



(6) 



In this expression (Equation 6), / is the power of the in- 
cident laser irradiation after the attenuation due to the dif- 
ferent optical elements in the light path, rj is the 
photothermal transduction efficiency (the parameter we 
want to calculate) that denotes a value for the efficacy of 
GNRs converting the incident light that interacts with 
them into thermal energy, and Ax is the optical density 
(also called absorbance) of the sample (colloidal disper- 
sion) at the irradiation wavelength. 

The outgoing heat flux can be considered linearly pro- 
portional to the thermal driving force, with a heat transfer 
coefficient, h, as proportionality constant: 



M(r-r,ef) 



(7) 



Therefore, the outgoing heat rate could be described 
using a lineal model with respect to the temperature, 
which results in the following equality when there is no 
incident laser light over the sample: 



dT 



Cp,i—= -hA{T-Tre{) 



(8) 



In the previous equations (Equations 7 and 8), r^f is 
the environment temperature and A is the irradiated area 
that the heat flux crosses toward the non-irradiated area. 
On the one hand, following this model, we can state that 
the part of the thermal cycle that defines the cooling of 
the sample exponentially depends on the time, and 
thereby, it is possible to determine the characteristic 
thermal time constant of the system by finding the 

Table 1 Values of the average time constant r,- and 
thermal capacitance for each tested volume 



Volume (^1) 




C„ (J/K) 


500 (/= 1) 


256.05 


13.06 


750 (/ = 2) 


295.15 


15.05 


1,000 (/ = 3) 


363.72 


18.55 
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Figure 5 Thermal capacitance values C,, (J/K) as a function of 
the sample volume (Vol). is the average squared error of the 
fitted line. 



0(W) = O.OlVol (|J|)+ 7.32 
R2 = 0.98 



exponential that adjusts the temperature curve. On the 
other hand, the heat transfer coefficient is inversely pro- 
portional to this time constant and could be defined as it 
is shown in the next expression: 



h- 



tA 



Q 

tA 



(9) 



Once we know the heat transfer coefficient, it can be 
used to calculate the amount of energy that the sample ac- 
cumulates or losses, from the temperature evolution. 

When the irradiated sample reaches the equilibrium 
temperature (T^^^), the incoming heat flux should be 
equal to the outgoing heat flux (toward the non-irradiated 
area of the sample), and therefore, r^jax is directly propor- 
tional to the optical power absorbed by GNRs: 



Ql + Qo- Qext = 0 

Q, + Qo = hA{T^,, 



*Ql + Qo = Qext 



(10) 



Therefore, we can solve the value of the photothermal 
transduction efficiency from the previous equations as 
follows: 



1 



hA{T^^^-T,^f)-QQ 



/(i-io- 



(11) 



Finally, if we rewrite the equation as a function of Q 
and rename T^^^ - T^gf = AT, we obtain the following 
expression, in which all the terms can be gotten from 
the calculation of the thermal parameters and the 
corresponding temperature curves for a fixed laser power 
irradiation: 



/(l-lO-^'i) 



(12) 



In this case, we have estimated the value of rj for the 
two kinds of nanorods described before (B-GNRs and 
PEG-GNRs). In both cases, the calculation was carried 
out from the temperature curves of three equal samples 
of water-dispersed GNRs with Ax = \ (which corresponds 
with a concentration of 36 |ig/ml) in order to obtain AT 



and from the temperature curves of three equal samples 
of deionized water to obtain Qq. All samples were irradi- 
ated with a laser power average of 2.0 W, and their vol- 
ume was 500 |il. 

Results and discussion 

Thermal parameters 

As described previously, we obtained three temperature 
curves of heating, stabilization, and cooling for each pro- 
posed case. Figure 3 shows schematically the shape of these 
curves and the parameters that we can get from them. 

We know that the thermal conductance is obtained 
from the data of power and temperature variation so that 
Cd = P / Ar„. Therefore, if we represent the value of P 
graphically as a function of AT^, it is possible to make a 
lineal fit in order to obtain the desired value of thermal 
conductance as shown in Figure 4. 

As it can be observed in Figure 4, the values of ther- 
mal conductance are pretty similar for the three consid- 
ered volumes. This behavior is consistent with the fact 
that the thermal conductance is an intensive magnitude, 
and therefore, it does not depend on the volume of the 
sample but on the global thermal properties of the 
considered system. 

Then, the thermal conductance of our system could be 
estimated from the average of the thermal conductances 
obtained for each volume: Qi (500 ^il) = 0.052 W/K, Q2 
(750 \A) = 0.052 W/K, and Q3 (1,000 ^1) = 0.048 W/K so 
that Q = 0.051 W/K. 

Then, Table 1 shows the average values of r, obtained 
for each tested volume and the associated values of 
thermal capacitance Q, (J/K), and Figure 5 represents 
graphically this evolution of the thermal capacitance as a 
function of the volume. 

In this case, we can observe that the values of thermal 
capacitance (expressed in J/K) have an increasing lineal 
evolution depending on the volume. This effect is consist- 
ent with the fact that the thermal capacitance is an exten- 
sive property, because their value depends on the amount 
of substance. Indeed, this parameter is defined as the ratio 




1500 2000 
Time (s) 

Figure 6 Temperature curves of heating, stabilization, and 
cooling (average values). Obtained for the irradiated samples of 
B-GNRs, PEG-GNRs, and deionized water (500 
'Aa = 1,Pla5er = 2.0 W). 
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Table 2 Parameters and photothermal transduction 
efficiency obtained from samples of B-BNRs and PEG- 
GNRs (500 Ml, A;^=^, Plaser = 2.0 W) 



/(W) 


GNRs 


Ar ro 


Qo(W) 


n (%) 


1.68 


B-GNRs 


20.5 


0.20 


56.0 




PEG-GNRs 


25.1 




714 



between the amount of thermal energy delivered to a 
body and the temperature change that the body experi- 
ments so that the greater is the volume, the lower is 
the experimented temperature change for a fixed 
amount of energy, and therefore, the greater is the 
value of the thermal capacitance. 

Photothermal transduction efficiency 

Once the thermal parameters have been calculated, to 
estimate the value of rj, it is needed to know the percent- 
age of the irradiation power (Plaser) incident on the 
sample (7) by taking into account the losses due to the 
light path through the optical elements between the out- 
put of the fiber and the sample. In this case, in the light 
path, there are a collimated lens and a 4-well plate. The 
lens has two faces and a percentage of losses of about 
3%, and the 4-well plate adds a percentage of losses near 
10% according to the manufacturers. Therefore, the total 
amount of losses is about 16%, and then, the incident 
power could be expressed as 7= (1 - 0.16) Plaser- 

Figure 6 shows the temperature curves obtained for 
the irradiated samples of B-GNRs, PEG-GNRs, and de- 
ionized water, Plaser = 2.0 W, and Table 2 shows the pa- 
rameters (average values) obtained from these 
temperature curves and the photothermal transduction 
efficiency that has been calculated from them according 
to Equation 12 (average values from three equal mea- 
sures in each case). 

Finally, Figure 7 shows the estimated values of tj for 
each one of the tested kinds of GNRs (B-GNRs and 
PEG-GNRs) presented as the average value (from the 
three temperature curves obtained for each case) ± the 
standard deviation. 

From the observation of the previous graph, we can 
affirm that under the established conditions of experi- 
mentation, the photothermal transduction efficiency of 
PEG-GNRs is about 15% to 16% higher than in B- 
GNRs. Thus, under the described working conditions, 
the use of PEG-GNRs allow for a more efficient heating 
of the sample, and therefore, the use of this kind of 
nanoparticles results in a more effective therapy. 

Conclusions 

The estimation of the thermal parameters and the 
photothermal transduction efficiency allows us to 
compare the thermal behavior of different kinds of 




B-GNRs PEG-GNRs 



Figure 7 Estimated values of r; for each one of the tested kinds 

of GNRs: B-GNRs and PEG-GNRs. Presented as the average value 

(from the three temperature curves obtained for each case) ±the 

standard deviation. 
^ 

nanoparticles under similar working conditions. More- 
over, it could be also useful to compare the thermal 
response of different kinds of nanoparticles under dif- 
ferent working conditions as, for example, concentra- 
tion, optical density, dispersion media, or sample 
holder. Therefore, the photothermal transduction effi- 
ciency is needed to determine the optimal conditions 
depending on each considered case. 

To summarize, we can say that, from a series of input 
data to the system, as the power of irradiation and the 
optical density of the used nanoparticles, it is possible to 
calculate the photothermal transduction efficiency of 
these particles using the thermal parameters of the sys- 
tem and the temperature variation of the samples. 
Therefore, it is possible to determine, for any kind of 
gold nanoparticles (or other noble metals) with their 
peak of absorption syntonized with the wavelength of ir- 
radiation, the percentage of the optical power that inter- 
acts (absorption + scattering) with the sample that really 
becomes in a temperature increasing. The higher the 
value of this parameter, the higher the efficiency of the 
designed optical hyperthermia treatment, and so, if we 
know the value of this parameter previously, we could 
select those nanoparticles that allow us to obtain better 
results in the designed therapy. 
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